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ABSTRACT 


We report the discovery of eight new ultra-faint dwarf galaxy candidates in the 
second year of optical imaging data from the Dark Energy Survey (DES). Six of these 
candidates are detected at high confidence, while two lower-confidence candidates are 
identified in regions of non-uniform survey coverage. The new stellar systems are found 
by three independent automated search techniques and are identified as overdensities of 
stars, consistent with the isochrone and luminosity function of an old and metal-poor 
simple stellar population. The new systems are faint {My > —4.7 mag) and span a 
range of physical sizes (17pc < ri /2 < 181 pc) and heliocentric distances (25kpc < 
Dq < 214kpc). All of the new systems have central surface brightnesses consistent 
with known ultra-faint dwarf galaxies (/U > 27.5 mag arcsec“^). Roughly half of the 
DES candidates are more distant, less luminous, and/or have lower surface brightnesses 
than previously known Milky Way satellite galaxies. Most of the candidates are found 
in the southern part of the DES footprint close to the Magellanic Clouds. We find 
that the DES data alone exclude {p < 10“^) a spatially isotropic distribution of Milky 
Way satellites and that the observed distribution can be well, though not uniquely, 
described by an association between several of the DES satellites and the Magellanic 
system. Our model predicts that the full sky may hold ~ 100 ultra-faint galaxies with 
physical properties comparable to the DES satellites and that 20-30% of these would 
be spatially associated with the Magellanic Clouds. 
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1. Introduction 

The population of Milky Way satellite galaxies includes the least luminous, least chemically en¬ 
riched, and most dark matter dominated galaxies in the known Universe (e.g., Simon & Geha 2007; 
Kirby et al. 2008). Although extreme systems from an observational perspective, low-luminosity 
dwarf spheroidal galaxies are likely to be the most common galaxy type by number. This duality 
places the emerging population of near-held galaxies in a unique position to test models of galaxy 
formation and the cold dark matter paradigm (e.g., Bullock et al. 2001), and has motivated studies 
ranging from the physical conditions at the time of reionization (e.g., Bullock et al. 2000; Benson 
et al. 2002) to the particle properties of dark matter (e.g., Ackermann et al. 2015; Geringer-Sameth 
et al. 2015a). 

Milky Way satellite galaxies are typically discovered as arcminute-scale statistical overdensities 
of individually resolved stars in wide-held optical imaging surveys (Willman 2010, and references 
therein). Prior to 2005, there were twelve known “classical” Milky Way satellite galaxies with 
absolute magnitudes in the range —8mag > My ^ —18mag (McGonnachie 2012). From 2005 
to 2014, hfteen ultra-faint satellite galaxies with My ^ —8 mag were identihed in Sloan Digital 
Sky Survey (SDSS; York et al. 2000) data through a combination of systematic searches (Willman 
et al. 2005a,b; Zucker et al. 2006a,b; Belokurov et al. 2006, 2007, 2008, 2009, 2010; Grillmair 2006, 
2009; Sakamoto & Hasegawa 2006; Irwin et al. 2007; Walsh et al. 2007) and dedicated spectroscopic 
follow-up observations (Kleyna et al. 2005; Munoz et al. 2006; Martin et al. 2007; Simon &: Geha 
2007; Geha et al. 2009; Koch et al. 2009; Belokurov et al. 2009; Walker et al. 2009; Garlin et al. 2009; 
Aden et al. 2009; Willman et al. 2011; Koposov et al. 2011; Simon et al. 2011; Kirby et al. 2013). 
Several outer halo star clusters and/or more compact stellar systems of uncertain classihcation were 
reported during roughly the same period (e.g., Koposov et al. 2007; Walsh et al. 2009; Belokurov 
et al. 2010; Munoz et al. 2012; Balbinot et al. 2013; Laevens et al. 2014; Belokurov et al. 2014; Kim 
& Jerjen 2015a; Kim et al. 2015b). Dwarf galaxies are distinguished from star clusters by having 
a dynamical mass that is substantially larger than the mass inferred from the luminous stellar 
population and/or a signihcant dispersion in stellar metallicities indicative of multiple generations 
of star formation and a gravitational potential deep enough to retain supernova ejecta (Willman & 
Strader 2012). 

Since the beginning of 2015, a combination of new optical imaging surveys and a reanalysis 
of SDSS data resulted in the detection of fourteen additional Milky Way companions, tentatively 
classihed as dwarf galaxies (Bechtol et al. 2015; Koposov et al. 2015a; Laevens et al. 2015a; Martin 
et al. 2015; Kim et al. 2015a; Kim & Jerjen 2015b; Laevens et al. 2015b). The new galaxy candidates 
have been examined for their (lack of) neutral gas (Westmeier et al. 2015), possible association with 
the Magellanic Glouds (Deason et al. 2015; Yozin &: Bekki 2015), spatial distribution within the 
Galactic halo (Pawlowski et al. 2015), and have become targets for indirect dark matter searches 
(Drlica-Wagner et al. 2015; Geringer-Sameth et al. 2015b; Hooper & Linden 2015). The provisional 
classification of these stellar systems as dwarf galaxies relies on their low surface brightnesses, large 
physical sizes, large ellipticities, and/or large heliocentric distances. However, some of these objects 
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lie in a region of size-luminosity space where the distinction between dwarf galaxy and globular 
cluster is ambiguous. Thus far, spectroscopic observations have confirmed that three of the recently 
discovered satellites, Reticulum II (Simon et al. 2015; Walker et al. 2015; Koposov et al. 2015b), 
Horologium I (Koposov et al. 2015b), and Hydra H (Kirby et al. 2015) possess the kinematic and 
chemical signatures of galaxies. 

The Dark Energy Survey (DES; Abbott et al. 2005), which began science operations in late 
2013, has already had a large impact on our census of Milky Way substructures. In the first year of 
DES data, nine new dwarf galaxy candidates were discovered in a region spanning less than 10% of 
the southern hemisphere (Bechtol et al. 2015; Koposov et al. 2015a; Kim &: Jerjen 2015b).^ Here, we 
present an extension of the analysis described in Bechtol et al. (2015) incorporating the second year 
of DES observations to expand the survey coverage from r\j 1800deg^ to ~5000deg^ (Eigure 1). In 
Section 2, we discuss the data reduction and catalog generation steps applied to the two-year DES 
data set and the resulting unique catalog of calibrated objects. We review the various algorithms 
applied to the DES data to search for ultra-faint galaxies in Section 3 and describe our candidate 
selection criteria in Section 4. The eight most significant stellar overdensities that are unassociated 
with previously known systems are reported in Section 5. If these objects are confirmed to be ultra¬ 
faint dwarf galaxies, they will be named after their resident constellations — Cetus H (DES J0117 
— 1725), Columba I (DES J0531—2801), Grus H (DES J2204—4626), Indus H (DES42038—4609),^ 
Reticulum HI (DES J0345—6026), Tucana HI (DES J2356—5935), Tucana IV (DES J0002—6051), 
and Tucana V (DES J2337—6316) — and named DES 2 through N if found to be globular clusters 
(e.g.. Tuque et al. 2015). We place the DES candidates in context with other known Local Group 
galaxies and the Magellanic system in Section 6, and we conclude in Section 7. 


2. Data Set 

DES is an ongoing optical imaging survey of ~ 5000 deg^ in the south Galactic cap using the 
Dark Energy Camera (DECam; Flaugher et al. 2015) on the 4-m Blanco Telescope at Cerro Tololo 
Inter-American Observatory (CTIO). The DECam focal plane comprises 62 2kx4k CCDs dedicated 
to science imaging and 12 2kx2k CCDs for guiding, focus, and alignment. The DECam field-of- 
view covers 3deg^ with a central pixel scale of Cf.'263. DES is scheduled for 525 nights distributed 
over five years, during which period each point in the survey footprint will be imaged ten times in 
each of the grizY bands (Abbott et al. 2005). 

The searches presented in Bechtol et al. (2015) and Tuque et al. (2015) utilized an object catalog 


^The system Kim 2/DES J2108.8—5109/Indus I was discovered independently in a separate data set by Kim et al. 
(2015b) slightly before the DES announcement. With deeper observations, Kim et al. concluded that this object is 
likely a star cluster, and we do not include it in our count of dwarf galaxy candidates. 

^To distinguish DES J2038—4609 from Kim 2/DES J2108.8—5109/Indus I, we suggest that the new system be 
designated Indus II if determined to be a dwarf galaxy. 
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Fig. 1.— Locations of the eight new dwarf galaxy candidates reported here (red triangles) along 
with nine previously reported dwarf galaxy candidates in the DES footprint (red circles; Bechtol 
et al. 2015; Koposov et al. 2015a; Kim & Jerjen 2015b), hve recently discovered dwarf galaxy 
candidates located outside the DES footprint (green diamonds; Laevens et al. 2015a; Martin et al. 
2015; Kim et al. 2015a; Laevens et al. 2015b), and twenty-seven Milky Way satellite galaxies known 
prior to 2015 (blue squares; McConnachie 2012). Systems that have been confirmed as satellite 
galaxies are individually labeled. The figure is shown in Galactic coordinates (Mollweide projection) 
with the coordinate grid marking the equatorial coordinate system (solid lines for the equator and 
zero meridian). The gray scale indicates the logarithmic density of stars with r < 22 from SDSS 
and DES. The two-year coverage of DES is ~5000deg^ and nearly fills the planned DES footprint 
(outlined in red). For comparison, the Pan-STARRS 1 Svr survey covers the region of sky with 
<^2000 > —30° (Laevens et al. 2015b). 
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generated from the coadded images in the DES year-one annual release (YlAl). To expedite the 
search for Milky Way satellites in year two, the present analysis uses data products derived from 
single-epoch imaging instead. This data set is referred to as the DES year-two quick release (Y2Q1), 
and its construction is summarized below. 


2.1. DES Year-Two Quick Release 

Observations: The Y2Q1 data set consists of 26,590 DECam exposures taken during the first 
two years of DES observing that pass DES survey quality cuts. Slightly less than half of the DES 
survey area was observed during the first season (Yl; 2013 August 15 to 2014 Eebruary 9) with 
typically two to four overlapping exposures, referred to as “tilings,” in each of the grizY filters. 
The second season (Y2; 2014 August 7 to 2015 February 15) covered much of the remaining survey 
area to a similar depth. Exposures taken in the griz bands are 90s, while Y-band exposures are 
45s. 

Image Reduction: The DES exposures were processed with the DES data management (DESDM) 
image detrending pipeline (Sevilla et al. 2011; Mohr et al. 2012; Desai et al. 2012, Gruendl et ah, 
in preparation). This pipeline corrects for cross-talk between CCD amplifier electronics, bias level 
variations, and pixel-to-pixel sensitivity variations (flat-fielding). Additional corrections are made 
for nonlinearity, fringing, pupil, and illumination. Both the Yl and Y2 exposures were reduced 
with the same image detrending pipeline used to process YlAl. 

Single-Epoch Catalog Generation: Astronomical source detection and photometry were per¬ 
formed on a per exposure basis using the PSFex and SExtractor routines (Bertin 2011; Berlin & 
Arnouts 1996) in an iterative process (Mohr et al. 2012, Gruendl et ah, in preparation). As part of 
this step, astrometric calibration was performed with SCAMP (Bertin 2006) by matching objects to 
the UCAC-4 catalog (Zacharias et al. 2013). The SExtractor source detection threshold was set 
at S/N > 10 for the Yl exposures, while for Y2 this threshold was lowered to S/N > 5 (the impact 
of this change on the stellar completeness is discussed in Section 2.2). During the catalog gener¬ 
ation process, we flagged problematic images (e.g., CCDs suffering from reflected light, imaging 
artifacts, point-spread function (PSF) mis-estimation, etc.) and excluded the affected objects from 
subsequent analyses. The resulting photometric catalogs were ingested into a high-performance 
relational database system. Throughout the rest of this paper, photometric fluxes and magnitudes 
refer to SExtractor output for the PSF model fit. 

Photometric Calibration: Photometric calibration was performed using the stellar locus regres¬ 
sion technique (SLR: Ivezic et al. 2004; MacDonald et al. 2004; High et al. 2009; Gilbank et al. 2011; 
Desai et al. 2012; Coupon et al. 2012; Kelly et al. 2014). Our reference stellar locus was empirically 
derived from the globally-calibrated DES YlAl stellar objects in the region of the YlAl footprint 
with the smallest E{B — V) value from the Schlegel et al. (SFD; 1998) interstellar extinction map. 
We performed a 1" match on all Yl and Y2 objects with S/N > 10 observed in r-band and at least 



one other band. We then applied a high-purity stellar selection based on the weighted average of the 
spreadjmodel quantity for the matched objects {\wavgspread.modeljr\ < 0.003; see below). The 
average zero point measured in YlAl, ZP^^izY = {30.0,30.3,30.2,29.9,28.0}, was assigned to each 
star as an initial estimate. Starting from this coarse calibration, we began an iterative procedure to 
fix the color uniformity across the survey footprint. We segmented the sky into equal-area pixels us¬ 
ing the HEALPix scheme (Gorski et al. 2005). For each ~0.2deg^ (resolution nside = 128) HEALPix 
pixel, we chose the DES exposure in each band with the largest coverage and ran a modified version 
of the Big MACS SLR code (Kelly et al. 2014)^ to calibrate each star from the reference exposure 
with respect to the empirical stellar locus. These stars became our initial calibration standards. 
We then adjusted the zero points of other CCDs so that the magnitudes of the matched detections 
agreed with the calibration set from the reference exposure. CCDs with fewer than 10 matched stars 
or with a large dispersion in the magnitude offsets of matched stars (azp > 0.1 mag) were flagged. 
For each calibration star, we computed the weighted-average magnitude in each band using these 
new CCD zero points; this weighted-average magnitude was used as the calibration standard for 
the next iteration of the SLR. In the first iteration, we assigned SLR zero points to the calibration 
stars based on the HEALPix pixel within which they reside. In subsequent iterations, we assigned 
SLR zero points to the calibration stars based on a bi-linear interpolation of their positions onto 
the HEALPix grid of SLR zero points. After the second iteration, the color zero points were stable 
at the 1 — 2 mmag level. The absolute calibration was set against the 2MASS J-band magnitude 
of matched stars (making use of the stellar locus in color-space), which were de-reddened using the 
SFD map with a reddening law of Aj = 0.709 x E{B — K)sfd from Schlafly & Finkbeiner (2011). 
The resulting calibrated DES magnitudes are thus already corrected for Galactic reddening by the 
SLR calibration. 

Unique Catalog Generation: There is significant imaging overlap within the DES footprint. 
For our final high-level object catalog, we identified unique objects by performing a (f!5 spatial 
match across all five bands. For each unique object, we calculated the magnitude in each hlter 
in two ways: (1) taking the photometry from the single-epoch detection in the exposure with the 
largest effective depth,^ and (2) calculating the weighted average {wavg) magnitude from multiple 
matched detections. During the catalog generation stage, we kept all objects detected in any two of 
the five filters (this selection is later restricted to the g- and r-bands for our dwarf galaxy search). 
The process of combining catalog-level photometry increases the photometric precision, but does 
not increase the detection depth. 

The resulting Y2Q1 catalog covers ~5100deg^ in any single band, ~5000deg^ in both g- and 
r-band, and ~4700deg^ in all five bands. The coverage of the Y2Q1 data set is shown in Figure 1. 
The Y2Q1 catalog possesses a median relative astrometric precision of ~ 40 mas per coordinate and 


^https://code.google.com/p/big-macs-calibrate/ 

"^The effective depth is calculated from the exposure time, sky brightness, seeing, and atmospheric transmission 
(Neilsen et al. 2015). 
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Fig. 2.— Performance of the morphological stellar classification based on \wavgspread-model-r\ 
< 0.003 + spreaderr-model-T was evaluated by comparing against deeper CFHTLenS imaging 
(Section 2.2). Left: Fraction of CFHTLenS-classified stars and galaxies that pass our DES stellar 
selection criteria. Right: Stellar purity evaluated in a typical high-Galactic latitude field centered 
at Galactic coordinates {l,b) = (176?3, —63?7). The shaded region represents the uncertainty 
associated with objects that are not confidently classified as either stars or galaxies with GFHTLenS. 
If a large fraction of these ambiguous objects are true stars, the inferred stellar purity for DES is 
near the upper bound of the indicated range. 

a median absolute astrometric uncertainty of ~ 140 mas per coordinate when compared against 
UGAG-4. The median absolute photometric calibration agrees within ~ 3% with the de-reddened 
global calibration solution of YlAl in the overlap region in the griz bands. When variations in the 
reddening law are allowed, this agreement improves to < 1%. The color uniformity of the catalog 
is ~ 1% over the survey footprint (tested independently with the red sequence of galaxies). The 
Y2Q1 catalog has a median point-source depth at a S/N = 10 of s' = 23.4, r = 23.2, i = 22.4, 
z = 22.1, and Y = 20.7. 


2.2. Stellar Object Sample 

We selected stars observed in both the g- and r-bands from the Y2Q1 unique object catalog 
following the common DES prescription based on the wavg spr ead-model and spr eaderr-model 
morphological quantities (e.g., Ghang et al. 2015; Koposov et al. 2015a). Specifically, our stellar 
sample consists of well-measured objects with \wavg spread-model-r\ < 0.003+spr eaderr-model-r, 
flags-{g,r} < 4, and magerr-psf-{g,r} < 1. By incorporating the uncertainty on spread-model, 
it is possible to maintain high stellar completeness at faint magnitudes while achieving high stellar 
purity at bright magnitudes. We sacrificed some performance in our stellar selection by using 
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the r-band spread-model which has a slightly worse PSF than the i-band; however, relaxing the 
requirement on i-band coverage increased the accessible survey area. 

The stellar completeness is the product of the object detection efficiency and the stellar clas¬ 
sification efficiency. As described in Section 2.1, different source detection thresholds were applied 
to the individual Y1 and Y2 exposures, and accordingly, the detection efficiency for faint objects in 
Y2Q1 varies across the survey footprint. We used a region of the DES Science Verification data set 
(SVAl), imaged to full survey depth (i.e., ~ 10 overlapping exposures in each of grizY), to estimate 
the object detection completeness using the individual Y1 exposures. For the Y2 exposures, we 
instead compared to publicly available catalogs from the Canada-France-Hawaii Telescope Lensing 
Survey (CFHTLenS; Erben et al. 2013). In both cases, the imaging for the reference catalog is 
~2mag deeper than the Y2Q1 catalog. These comparisons show that the Y2Q1 object catalog is 
> 90% complete at about 0.5 mag brighter than the S/N = 10 detection limit. 

We estimated the stellar classification efficiency by comparing to CFHTLenS and by us¬ 
ing statistical estimates from the Y2Q1 data alone. For the CFHTLenS matching, we identi¬ 
fied probable stars in the CFHTLenS catalog using a combination of the classstar output from 
SExtractor and the fitclass output from the lensfit code (Heymans et al. 2012). Specifically, 
our CFHTLenS stellar selection was {fitclass = 1) OR {cl ass star > 0.98), and our galaxy selec¬ 
tion was {fitclass = 0) OR {cl ass star < 0.2). Comparing against the DES stellar classification 
based on spreadjmodel, we found that the Y2Q1 stellar classification efficiency exceeds 90% for 
r < 23 mag, and falls to ~80% by r = 24 (Figure 2). We independently estimated the stellar 
classification efficiency by creating a test sample of Y2Q1 objects with high stellar purity using the 
color-based selection r — i > 1.7. We then applied the spread-model morphological stellar selection 
to evaluate the stellar completeness for the test sample. The stellar selection efficiency determined 
by the color-based selection test is in good agreement with that from the CFHTLenS comparison. 

For a typical field at high Galactic latitude, the stellar purity is ~ 50% at r = 22.5 mag, 
and falls to ~ 10% at the faint magnitude limit (Figure 2). In other words, the majority of faint 
objects in the “stellar sample” are misclassified galaxies rather than foreground halo stars.® This 
transition can be seen in the color-magnitude distribution diagnostic plots for individual candidates 
(Section 6.1). 


3. Search Methods 

In this section we describe three complementary search strategies used to identify overdensities 
of individually resolved stars matching the old and metal-poor populations characteristic of ultra¬ 
faint galaxies. Each technique used a different method to estimate the spatially varying field 


®This conclusion differs from that of Koposov et al. (2015a), who state that low levels of contamination are 
observed to the magnitude limit of DES. 
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density and to select potential satellite member stars. The output of each search method was a set 
of seed coordinates, ( 02000 ; <^ 2000 ; M), corresponding to statistical excesses relative to the local 

field density. These seeds were cross-matched between methods and examined in greater detail 
(Section 5). 

When searching for new stellar systems, each detection algorithm applied a deep magnitude 
threshold {g < 23.5 mag or g < 24mag depending on the search technique). At these faint mag¬ 
nitudes, our stellar sample is far from complete (Section 2.2). However, real stellar systems are 
dominated by faint stars, and the signal-to-noise ratio for satellite discovery can increase by in¬ 
cluding objects that are fainter than the stellar completeness limit. After candidate systems were 
identified, we applied a brighter magnitude threshold {g < 23 mag) to avoid biasing the fit results 
during object characterization (Section 5). 

One potential concern when extending our search to very faint magnitudes is the increased 
contamination from misclassified galaxies (Section 2.2). However, as long as the contaminating 
objects have a similar statistical distribution in the target region and the surrounding area, each 
technique naturally incorporates the increased object density into the background estimate. There 
are situations where the assumption of uniformity does not apply (e.g., clusters of galaxies, varia¬ 
tions in survey depth, gaps in coverage, etc.); therefore, when assembling our candidate list, we also 
examined the distribution of probable galaxies as well as stars with colors far from the expected 
isochrone (Section 5). 

We briefly describe our three search strategies below. 


3.1. Stellar Density Maps 

Our most straightforward and model-independent search technique used a simple isochrone 
filter to facilitate visual inspection of the stellar density field. The specific algorithm was a variant 
of the general approach described in Section 3.1 of Bechtol et al. (2015) and follows from the 
methods of Koposov et al. (2008), Walsh et al. (2009), and Kim h Jerjen (2015a). 

We first increased the contrast of putative satellite galaxies by selecting stars that were con¬ 
sistent with the isochrone of an old (r = 12 Gyr) and metal-poor {Z = 0.0001) stellar population 
(Dotter et al. 2008). Specifically, we selected stars within 0.1 mag of the isochrone locus in color- 
magnitude space, enlarged by the photometric measurement uncertainty of each star added in 
quadrature. After isochrone filtering, we smoothed the stellar density field with a 2' Gaussian 
kernel. For each density peak in the smoothed map, we scanned over a range of radii (1^ to 18') 
and computed the Poisson significance of finding the observed number of stars within each circle 
centered on the density peak given the local field density. This process was repeated for stellar 
populations at a range of distance moduli (16 mag < m — M < 24 mag) and diagnostic plots were 
automatically generated for each unique density peak. Our nominal magnitude threshold for this 
search was 5 ^ < 24 mag; however, in regions on the periphery of the survey with large variations in 
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coverage we used a brighter threshold of 5 < 23 mag. 


3.2. Matched-filter Search 

Our second search strategy utilized a matched-filter algorithm (Luque et al. 2015; Balbinot 
et al. 2011, 2013). We began by binning objects in our stellar sample with 17mag < g < 24mag into 
1' X V spatial pixels and g versus g — r color-magnitude bins of 0.01 mag x 0.05 mag. We then created 
a grid of possible simple stellar populations modeled with the isochrones of Bressan et al. (2012) 
and populated according to a Chabrier (2001) initial mass function to use as signal templates. The 
model grid spanned a range of distances (lOkpc < D < 200kpc), ages (9 < log^g < 10-2), and 
metallicities (Z = {0.0002,0.001,0.007}). The color-magnitude distribution of the field population 
was empirically determined over larger 10° x 10° regions. For each isochrone in the grid, we then 
fit the normalization of the simple stellar population in each spatial pixel to create a corresponding 
density map for the full Y2Q1 area. 

The set of pixelized density maps do not assume a spatial model for the stellar system, and can 
be used to search for a variety of stellar substructures (e.g., streams and diffuse clouds). To search 
specifically for dwarf galaxies, we convolved the density maps with Gaussian kernels having widths 
ranging from 0^ (no convolution) to 9' and used SExtractor to identify compact stellar structures 
in those maps. This convolution and search was performed on each map for our grid of isochrones. 
The resulting seeds were ranked according to their statistical significance and by the number of 
maps in which they appeared, i.e., the number of isochrones for which an excess was observed. 


3.3. Likelihood-based Search 

Our third search technique was a maximum likelihood-based algorithm that simultaneously 
used the full spatial and color-magnitude distribution of stars, as well as photometric uncertainties 
and information about the survey coverage. The likelihood was constructed from the product 
of Poisson probabilities for detecting each individual star given its measured properties and the 
parameters that describe the field population and putative satellite (Section 3.2 of Bechtol et al. 
2015). When searching for new stellar systems, we used a spherically symmetric Plummer model 
as the spatial kernel, and a composite isochrone model consisting of four isochrones of different 
ages, r = {12 Gyr, 13.5 Gyrj, and metallicities, Z = {0.0001, 0.0002}, to bracket a range of possible 
stellar populations (Bressan et al. 2012). We scanned the Y2Q1 data testing for the presence 
of a satellite galaxy at each location on a multi-dimensional grid of sky position (0l7 resolution; 
HEALPix nside = 4096), distance modulus (16mag < m — M < 23mag), and kernel half-light 
radius (rh = {0?03, 0?1, 0?3}). For the likelihood search, we restricted our magnitude range to 
16 mag < g < 23.5 mag as a compromise between survey depth and stellar completeness. The 
statistical significance of a putative galaxy at each grid point was expressed as a test statistic (TS) 
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based on the likelihood ratio between a hypothesis that includes a satellite galaxy versus a field- 
only hypothesis (Equation 4 of Bechtol et al. 2015). As a part of the model fitting, the satellite 
membership probability is computed for every star in the region of interest. These photometric 
membership probabilities can be used to visualize which stars contribute to the candidate detection 
significance (Section 5) and have previously been used to select targets for spectroscopic follow¬ 
up (Simon et al. 2015; Walker et al. 2015). 


4. Candidate Selection 

The search methods described in Section 3 each produced a set of seed locations ( 02000 ; <^ 2000 ; 
m — M) for overdensities in the stellar field and a significance associated with each seed. Before 
performing a computationally-intensive multi-parameter likelihood fit to characterize each seed, 
we applied a set of simple selection cuts to remove seeds that are unlikely to be new Milky Way 
satellites. We set detection significance thresholds at > 5.5 ct for the stellar density map search 
and TS > 45 (~ 6cr) for the maximum-likelihood search. For the matched-filter method, the ten 
highest-ranked seeds were selected from each 10° x 10° search region. After thresholding, the 
union of all three search techniques yielded several hundred unique seeds. Many of these seeds 
were attributed to steep gradients in the stellar density field, numerical effects near the survey 
boundaries, and imaging artifacts. We also compared against catalogs of other astrophysical objects 
that produce false positives, such as large nearby galaxies (Nilson 1973; Corwin 2004) and galaxy 
clusters (Rykoff et al. 2014). These objects were removed from our list of candidates and we did 
not pursue investigation at lower significance thresholds due to an increased false positive rate. 

The resulting list of seeds was matched against catalogs of known Local Group galaxies (Mc- 
Connachie 2012) and star clusters (Kharchenko et al. 2013; Harris 1996, 2010 edition). Some of 
the most prominent stellar systems in our list of seeds were the Fornax, Phoenix, Sculptor,® and 
Tucana dwarf galaxies, and the globular clusters AM 1, Eridanus, Reticulum, Whiting 1, NGC 288, 
NGC1261, NGG1851, NGC 1904, and NGC 7089. Additionally, seed locations were compared 
against other stellar overdensities recently reported in the DES footprint (Kim et al. 2015b; Bech¬ 
tol et al. 2015; Koposov et al. 2015a; Kim & Jerjen 2015b; Luque et al. 2015) and all eleven objects 
were recovered with high significance. The locations of previously known stellar systems in the 
DES footprint are shown in Figure 3. 

As an additional test of the physical nature of the ~ 30 highest significance unassociated seeds, 
we explored the posterior probability density of the likelihood with an affine invariant Markov Chain 
Monte Carlo (MCMC) ensemble sampler (MCMC; Foreman-Mackey et al. 2013).'^ This step was 
meant as a preliminary investigation of whether the posterior distribution was well constrained by 


°The center of Sculptor lies in a hole in Y2Q1 coverage, though the outskirts are still detected at high significance, 
^emcee v2.1.0: http://dan.iel.fm/emcee/ 
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Fig. 3.— Cartesian projection of the density of stars observed in both g- and r-bands with g < 23 
and g — r < 1 over the DES Y2Q1 footprint (~ 5000 deg^). Globular clusters are marked with 
symbols (Harris 1996, 2010 edition), two faint outer halo clusters are marked with “x” symbols 
(Kim et al. 2015b; Luque et al. 2015), Local Group galaxies known prior to DES are marked with 
blue squares (McConnachie 2012), dwarf galaxy candidates discovered in Y1 DES data are marked 
with red outlined circles, while the Y2 stellar systems are marked with red triangles. The periphery 
of the LMC can be seen in the southeast corner of the footprint, while the Galactic stellar disk can 
be seen on the eastern and western edges. 
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the data. We fixed the distance modulus of each object at the value that maximized the likelihood in 
the grid search (Section 3.3) and sampled the morphological parameters (center position, extension, 
ellipticity, and position angle) along with the richness of the putative stellar system. We imposed 
a Jeffreys’ prior on the extension {V{ah) oc 1/ah) and flat priors on all other parameters. The 
posterior probability was sampled with 2.5 x 10® steps and we rejected seeds that lacked a well- 
defined maximum in the posterior. 

The process described above resulted in the selection of eight candidate stellar systems. Six of 
these candidates, DES J2204—4626 (Gru I), DES J2356—5935 (Tuc III), DES J0531—2801 (Col I), 
DES J0002—6051 (Tuc IV), DES J0345—6026 (Ret III), and DES J2337—6316 (Tuc V), are high- 
confidence detections in clean regions of the survey. The additional two candidates, DES J2038 
—4609 (Ind II) and DES J0117—1725 (Get II), are lower confidence and reside in more complicated 
regions. In Table 1 we report the coordinates and detection significances of each of these objects. 


5. Candidate Characterization 

In this section, we describe the iterative procedure used to characterize each of the eight candi¬ 
date stellar systems. When fitting the new candidates, we applied a brighter magnitude threshold, 
g < 23 mag, to mitigate the impact of stellar incompleteness. The results of our characterization 
are shown in Table 1 and Table 2. 

We began by simultaneously fitting the morphological parameters and distance modulus of 
each candidate following the procedure described in Section 4. Best-fit values were derived from 
the marginalized posterior distribution and the morphological parameters were temporarily hxed at 
these values. We then ran a MGMC chain simultaneously sampling the distance, age, metallicity, 
and richness, assuming flat priors on each parameter. Signihcant correlations between these param¬ 
eters were found, and in some cases the age and metallicity were poorly constrained (see below). To 
assess the error contribution from intrinsic uncertainty in the isochrone, we resample the posterior 
distribution using the Dartmouth isochrone family (Dotter et al. 2008). In general, the best-fit 
distance moduli agree to within 0.1 mag when the data were fit by these two isochrone classes. To 
determine the uncertainty on the distance modulus quoted in Table 1 we started with the highest 
posterior density interval from the Bressan et al. (2012) isochrone fit. We then calculated a sys¬ 
tematic component from the difference in the best-fit distance modulus derived with Bressan et al. 
(2012) and Dotter et al. (2008) isochrones. An additional ±0.03 mag uncertainty was added to 
account for uncertainty in the photometric calibration. The age and metallicity values were taken 
from the peak of the marginalized posterior density from the Bressan et al. (2012) isochrone fit. 

Eor the distant systems (i.e., DES J0531—2801, DES J2038—4609, and DES J0345—6026), the 
main-sequence turnoff (MSTO) is fainter than the g < 23 mag limit imposed for our fit. In these 
systems, the data provide weaker constraints on the age and metallicity. We therefore fit the 
distance modulus using a composite of four Bressan et al. (2012) isochrones with hxed ages, r = 
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{12 Gyr, 13.5 Gyr}, and metallicities, Z = {0.0001,0.0002}. We followed the procedure above to 
incorporate systematic uncertainties from the choice of theoretical isochrone family and from the 
calibration.® 

Fixing the distance, age, and metallicity at the values derived in the previous step, we then 
re-fit the morphological parameters. The best-fit values of the morphological parameters and 
their highest posterior density intervals are listed in Table 1. The majority of objects show no 
significant evidence for ellipticity, which is confirmed by calculating the Bayes factor (BF) via the 
Savage-Dickey density ratio (Dickey 1971; Trotta 2007). For objects with evidence for asphericity 
{BF > 3), we report both the elliptical Plummer radius, ah, and the position angle, (p. For all 
objects, we report the azimuthally averaged half-light radius, Vh- In the region of DES J2356—5935 
there is a linear structure visible in the filtered stellar density map, consistent with a set of tidal 
tails (Section 6.1). We require the spatial profile of DES J2356—5935 to be azimuthally symmetric 
in the MCMC fit and analyze this linear feature separately. 

The best-fit parameters derived from our iterative MGMC analysis along with several addi¬ 
tional derived properties are reported in Table 1 and Table 2. The physical size of each object was 
calculated by propagating the uncertainty from the distance and azimuthally averaged half-light 
radius. Absolute magnitudes were calculated according to the prescription of Martin et al. (2008) 
and do not include the distance uncertainty. 

The two lower confidence candidates are located in complicated regions of the survey. DES J0117 
— 1725 is located in a region of sparse coverage and is close to both a CCD chip gap and a larger 
hole in the survey. DES J2038—4609 is located in a region of non-uniform depth at the interface of 
Y1 and Y2 observations. Although we have attempted to account for these issues in the likelihood 
fit, we caution that the parameters derived for these systems may be less secure than for the other 
candidates. 

As a final note, the iterative MCMC fitting procedure described above was also applied to the 
satellite galaxy candidate Crus I (Koposov et al. 2015a). Crus I was identified in YlAl imaging 
data, but was located close to a CCD chip gap in a region with sparse coverage. We reanalyzed 
Crus I with the additional Y2 exposures and found that its structural parameters are consistent 
within uncertainties with those reported by Koposov et al. (2015a). 


6. Discussion 


Spectroscopic observations are needed to definitively classify the newly discovered stellar sys¬ 
tems as either star clusters or ultra-faint dwarf galaxies based upon their stellar dynamics and/or 


®One notable exception is the case of DES J2038—4609, where significant constraining power conies from a set of 
horizontal branch (HB) stars. For DES J2038—4609 we estimated the distance uncertainty from the Bressan et al. 
(2012) isochrone analysis alone. 
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Table 1. Detection of ultra-faint galaxy candidates. 


Name 

c ^2000 

(deg) 

<52000 

(deg) 

m — AT 

f'h 

C) 

C) 

€ 

(f> 

Map Sig. 

TS Scan 

EPi 

DES J2204-4626 (GruII) 

331.02 

-46.44 

18.62 ± 0.21 

0 q+0.9 
'^■^-0.5 


< 0.2 


15.7 

369 

161 

DES J2356-5935 (TucIII)^ 

359.15 

-59.60 

17.01 ± 0.16 

0 q+ 0-8 




11.1 

390 

168 

DES J0531-2801 (Coll)'^ 

82.86 

-28.03 

21.30 ± 0.22 

1 9+0-5 

^■^-0.4 


< 0.2 


10.5 

71 

33 

DES J0002-6051 (TucIV) 

0.73 

-60.85 

18.41 ± 0.19 

9 1 + ^-'^ 


6 6 
+ 1 

6 

11 ±9 

8.7 

287 

134 

DES J0345-6026 (Retlll)*^ 

56.36 

-60.45 

19.81 ± 0.31 

2 4+0-0 
■ - 0.8 


< 0.4 


8.1 

56 

22 

DES J2337-6316 (TucV) 

354.35 

-63.27 

18.71 ± 0.34 

1 0+O-3 

^-*-^-0.3 


0 7+01 
‘ - 0.2 

30 ± 5 

8.0 

129 

24 

DES J2038-4609 (Indll)'^ 

309.72 

-46.16 

21.65 ± 0.16 

2 9+^-^ 


< 0.4 


6.0 

32 

22 

DES J0117-1725 (Get II) 

19.47 

-17.42 

17.38 ± 0.19 

1 9+1-0 

^■^-0.5 


< 0.4 


5.5 

53 

21 


Note. — Characteristics of satellite galaxy candidates discovered in DES Y2 data. Best-fit parameters from the maximum-likelihood analysis 
assume a Bressan et al. (2012) isochrone. Uncertainties come from the highest density interval containing 68 % of the posterior distribution. 
The uncertainty on the distance modulus (m — AT) also includes systematic uncertainties coming from the choice of theoretical isochrone and 
photometric calibration (Section 5). The azimuthally averaged half-light radius (r/^) is quoted for all candidates. For systems with evidence for 
asphericity (Bayes’ factor > 3), we quote the ellipticity (e), the position angle (</>), and the length of the semi-major axis of the ellipse containing 
half of the light (a;^ = /\/l — e). Upper limits on the ellipticity are quoted for other candidates at 84% confidence. “Map Sig.” refers to 
detection significance of the candidate from the stellar density map search method (Section 3.1). “TS Scan” refers to the significance (Equation 
4 in Bechtol et al. 2015) from the likelihood scan (Section 3.3). Ep^ is the estimated number of satellite member stars in the DES Y2Q1 catalog 
with g < 23 mag. 

^Fit with a spherically symmetric Plummer profile due to the possible presence of tidal tails (Section 6.1). 

’^Fit with a composite isochrone: r = {12 Gyr, 13.5 Gyr}, Z = {0.0001, 0.0002} (Section 5). 


Table 2. Derived properties of ultra-faint galaxy candidates. 


Name T h Distance AT* AT^ '^1/2 ^ 




(deg) 

(deg) 

(kpc) 

(io^Mq) 

(mag) 

(pe) 

(Gyr) 


DES J2204-4626 

DES J2356-5935 

DES J0531-2801 

DES J0002-6051 

DES J0345-6026 

DES J2337-6316 

(GruII) 

(TucIII)^ 

(GolI)O 

(TucIV) 

(Retlll)O 

(TucV) 

351.15 

315.38 

231.62 

313.29 

273.88 

316.31 

-51.94 

-56.19 

-28.88 

-55.29 

-45.65 

-51.89 

53 ± 5 

25 ± 2 

182 ± 18 

48 ± 4 

92 ± 13 

55 ± 9 

o. 8 «;i 

0 2 + ^-0 

9 9+0.4 
^•^-0,3 

2 q+ 0.6 

^•0-0.7 

0 5+0-1 

-3.9 ± 0.22 

-2.4 ± 0.42 

-4.5 ± 0.17 

-3.5 ± 0.28 

-3.3 ± 0.29 

-1.6 ± 0.49 

93 ± 14 

44 + 6 

103 ± 25 

127 ± 24 

64 ± 24 

17 + 6 

12.5 

10.9 

11.6 

10.9 

0.0002 

0.0001 

0.0001 

0.0003 

DES J2038-4609 (Indll)^ 
DES J0117-1725 (Get II) 

353.99 

156.48 

-37.40 

-78.53 

214 ± 16 

30 ± 3 

4 9+1-® 

0 1+0.04 
^•-^-0.04 

-4.3 ± 0.19 

0.0 ± 0.68 

181 + 67 

17 + 7 

10.9 

0.0002 


Note. — Derived properties of the DES Y2 satellite candidates. Stellar masses (AT*) are computed for a Chabrier 
(2001) initial mass function. The absolute visual magnitude is derived via the procedure of Martin et al. (2008) using the 
transformation equations from Bechtol et al. (2015). The uncertainty on the azimuthally-averaged physical half-light radius 
(ri/ 2 ) includes the uncertainty in the projected half-light radius and distance. Age (r) and metallicity {Z) values come 
from the peak of the posterior distribution. 

^Fit with a spherically symmetric Plummer profile due to the possible presence of tidal tails (Section 6.1). 

'^Fit with a composite isochrone: r = {12 Gyr, 13.5 Gyr}, Z = {0.0001, 0.0002} (Section 5). 
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Fig. 4.— Local Group galaxies (McConnachie 2012) and globular clusters (Harris 1996, 2010 
edition) occupy distinct regions in the plane of physical half-light radius (azimuthally averaged) 
and absolute magnitude. The majority of DES satellite candidates (red triangles and circles) are 
more consistent with the locus of Local Group galaxies (blue squares) than with the population 
of Galactic globular clusters (black “-I-” ). Other recently reported dwarf galaxy candidates (green 
diamonds) include Hydra II (Martin et al. 2015), Triangulum H (Laevens et al. 2015a), Pegasus HI 
(Kim et al. 2015a), Draco II and Sagittarius H (Laevens et al. 2015b). Several outer halo star 
clusters and systems of ambiguous classification are indicated with “x” symbols: Koposov 1 and 
Koposov 2 (Koposov et al. 2007; Paust et al. 2014), Segue 3 (Belokurov et al. 2010; Fadely et al. 
2011; Ortolan! et al. 2013), Munoz 1 (Munoz et al. 2012), Balbinot 1 (Balbinot et al. 2013), 
Laevens 1 (Laevens et al. 2014; Belokurov et al. 2014), Laevens 3 (Laevens et al. 2015b), Kim 1 
and Kim 2 (Kim h Jerjen 2015a; Kim et al. 2015b), and DES 1 (Luque et al. 2015). Dashed lines 
indicate contours of constant surface brightness at ^ = {25, 27.5, 30} mag arcsec“'. 
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metallicity dispersions. However, the distances, physical sizes, and luminosities derived from DES 
photometry (Table 2) already suggest a galactic classification for many of the new candidates. In 
Figure 4, we show the distribution of physical half-light radius versus absolute magnitude for Milky 
Way globular clusters (Harris 1996, 2010 edition) and Local Group galaxies (McConnachie 2012). 
We find that all of the systems discovered in Y2Q1 fall along the established locus for nearby galax¬ 
ies. Several of the new systems possess lower surface brightnesses than any conhrmed ultra-faint 
galaxy,® supporting earlier conclusions that the threshold in surface brightness was an observational 
selection effect (Koposov et al. 2008; Walsh et al. 2009). Finally, globular clusters generally have 
small ellipticities, e < 0.2 (van den Bergh 2008), while dwarf galaxies are commonly more elliptical 
(Martin et al. 2008). However, the ellipticities of many of the new objects are not well constrained 
by the DES data and thus do not provide a strong indicator for the nature of these objects. 

The spatial distribution of the new candidates within the DES footprint is suggestive of an 
association with the Magellanic Clouds (Section 6.3). When discussing this scenario, the coordinates 
of the Large Magellanic Cloud (LMC) are taken to be (a 2000 ) ^ 2000 ) Dq) = (80?89, —69?76,49.89 kpc) 
(distance from de Grijs et al. 2014) and the coordinates of the Small Magellanic Cloud (SMC) are 
taken to be (a 2000 ) <^20005 Dq) = (13?19, —72?83, 61.94kpc) (distance from de Grijs & Bono 2015).^® 
While the new DES candidates reside in a region close to the Magellanic Stream, no candidate is 
found to be coincident with the main filament of the Stream. A similar observation was made by 
Koposov et al. (2015a) for the DES stellar systems discovered in Yl. 

We also investigate a potential group of satellites in the constellation Tucana, consisting of 
DES J0002—6051 and DES J2337—6316, and the object DES J2251.2—5836 (Tuc H) discovered 
in Yl data (Bechtol et al. 2015; Koposov et al. 2015a). The centroid of the Tucana group is 
at (0:2000) <52000) .Dq) = (351?90,—61?06, 53.63 kpc), and the separation of each member from the 
centroid is < 7 kpc (Table 3). This grouping of objects is projected onto a region of H i gas that 
is likely a secondary filament or high-velocity cloud associated with the head of the Magellanic 
Stream (Putman et al. 2003; McClure-Griffiths et al. 2009). 


6.1. Review of Individual Candidates 

• DES J2204—4626 (Crus H, Figure 6 ): As the most significant new candidate, DESJ2204 
—4626 has ~ 161 probable member stars with g < 23 mag in the DES imaging. The large 
physical size of this system (93 pc) allows it to be tentatively classihed as an ultra-faint dwarf 
galaxy. A clear red giant branch (RGB) and MSTO are seen in the color-magnitude diagram 
of DES J2204—4626. Several likely blue horizontal branch (BHB) members are seen at <7 ~ 


®A potential exception to this statement is the large, faint, and perhaps tidally disrupted object Bootes III (Carlin 
et al. 2009). 

^°Sky coordinates taken from NED: https://ned.ipac.caltech.edu. 
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Fig. 5.— DES satellite galaxy candidates (red circles and triangles) have comparable luminosities to 
known Local Group galaxies (blue squares; McConnachie 2012) and other newly discovered galaxy 
candidates (green diamonds), but are typically found at larger distances. Dashed lines indicate the 
heliocentric distances of the LMC and SMC. Globular clusters (black “+” symbols, Harris 1996) 
and outer halo star clusters (“x” symbols) have been detected at comparable luminosities and 
distances due to their more compact nature. The constituents of each source class can be found in 
Figure 4. 
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19 mag. There is a small gap in survey coverage ~ 0?5 from the centroid of DES J2204—4626, 
but this is accounted for in the maximum-likelihood analysis and should not significantly 
affect the characterization of this rich satellite. 

DES J2204—4626 is nearly equidistant from the LMC (46kpc) and from the Galactic center 
(49 kpc). It is slightly closer to the SMC (33kpc) and intriguingly close to the the group of 
stellar systems in Tucana (~ 18kpc). While it is unlikely that DES J2204—4626 is currently 
interacting with any of the other known satellites, we cannot preclude the possibility of a past 
association. 

• DES J2356—5935 (Tucana III, Figure 7): DES J2356—5935 is another highly signihcant 
stellar system that presents a clearly defined main sequence and several RGB stars. The 
physical size and luminosity of DES J2356—5935 (ri /2 = 44pc. My = —2.4mag) are compa¬ 
rable to that of the recently discovered dwarf galaxy Reticulum II. At a distance of 25 kpc, 
DES J2356—5935 would be among the nearest ultra-faint satellite galaxies known, along with 
Segue 1 (23kpc), Reticulum II (32 kpc), and Ursa Major II (32kpc). The relative abundance 
of bright stars in DES J2356—5935 should allow for an accurate spectroscopic determination 
of both its velocity dispersion and metallicity. DES J2356—5935 is reasonably close to the 
LMG and SMC (32 kpc and 38kpc, respectively), and measurements of its systemic velocity 
will help to elucidate a physical association. 

As mentioned in Section 5, there is an additional linear feature in the stellar density field 
around DES J2356—5935 that strongly influences the fitted ellipticity of this galaxy candi¬ 
date. In Figure 8, we show the 6° x 6° region centered on DES J2356—5935 which contains 
this linear feature. Selecting stars with the same isochrone filter used to increase the con¬ 
trast of DES J2356—5935 relative to the held population, we found a faint stellar overdensity 
extending ~2° on both sides of DES J2356—5935. This feature has a FWHM of ~0?3 (i.e., 
projected dimensions of 1.7 kpc x 130 pc at a distance of 25 kpc) and is oriented ~85° east 
of north. The length of the feature and the smooth density held observed with an inverted 
isochrone hlter support the conclusion that this is a genuine stellar structure associated with 
DES J2356-5935. 

One interpretation is that DES J2356—5935 is in the process of being tidally stripped by 
the gravitational potential of the Milky Way. The relatively short projected length of the 
putative tails may imply that DES J2356—5935 is far from its orbital pericenter (where tidal 
effects are strongest) and/or that stripping began recently. Despite the possible presence of 
tidal tails, the main body of DES J2356—5935 appears relatively round. While the proximity 
of DES J2356—5935 might make it an important object for indirect dark matter searches, 
evidence of tidal stripping would suggest that a large fraction of its outer dark matter halo 
has been removed. However, the mass within the stellar core (and therefore the dark matter 
content in the central region) can still likely be determined accurately (Oh et al. 1995; Munoz 
et al. 2008; Peharrubia et al. 2008). 


DES J0531—2801 (Columba I, Figure 9): DES J0531—2801 is the second most distant of the 
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DES J2204-4626 (Grus II); (aaooo, <^ 2000 , m-M) = (33E02, -46?44, 18.62) 




g-r 


Fig. 6.— Stellar density and color-magnitude diagrams for DES J2204—4626 (Grus II). Top left: 
Spatial distribution of stellar objects with g < 24mag that pass the isochrone filter (see text). The 
held of view is 2° x 2° centered on the candidate and the stellar distribution has been smoothed 
by a Gaussian kernel with standard deviation 1'2. Top center: Radial distribution of objects 
with g — r < Imag and g < 24mag: stars passing the isochrone hlter (red), stars excluded from 
the isochrone hlter (gray), and galaxies passing the isochrone hlter (black). Top right: Spatial 
distribution of stars with high membership probabilities within a 0?5 x 0?5 held of view. Small 
gray points indicate stars with membership probability less than 5%. Bottom left: Same as top left 
panel, but for probable galaxies that pass the isochrone hlter. Bottom center: The color-magnitude 
distribution of stars within 0?1 of the centroid are indicated with individual points. The density of 
the held within an annulus from radii of 0?5 to 1?0 is represented by the background two-dimensional 
histogram in grayscale. The blue curve shows the best-ht isochrone as described in Table 1 and 
Table 2. Bottom right: Color-magnitude distribution of high membership probability stars. 
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Table 3. Physical separations of ultra-faint galaxy candidates. 


Name 


DgC 

^LMC 

-Dsmc 

D'T'iir 


(kpc) 

(kpc) 

(kpc) 

(kpc) 

(kpc) 

DES J2204-4626 (GruII) 

53 

49 

46 

33 

18 

DES J2356-5935 (TucIII) 

25 

23 

32 

38 

28 

DES J0531-2801 (Coll) 

182 

186 

149 

157 

168 

DES J0002-6051 (Tuc IV) 

48 

46 

27 

18 

7 

DES J0345-6026 (Ret III) 

92 

92 

45 

40 

53 

DES J2337-6316 (TucV) 

55 

52 

29 

14 

3 

DES J2038-4609 (Indll) 

214 

208 

193 

170 

169 

DES J0117-1725 (Cet II) 

30 

32 

46 

51 

40 

DES J2251.2-5836 (Tuc II) 

58 

54 

37 

20 

7 


Note. — Three-dimensional physical separation between DES satellite candi¬ 
dates and the LMC and SMC. Also listed is the heliocentric distance, the Galac- 
tocentric distance, and the separation to the centroid of the proposed Tucana 
group comprising DES J2251.2-5836, DES J0002-6051, and DES J2337-6316. 


DES J2356-5935 (Tucana III); (azooo, ^ 2000 , m - M) = (359?15, -59?60, 17.01) 



Fig. 7. 


Analogous to Figure 6 but for DES J2356—5935 (Tucana III). 
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(a 2000 . ^ 2000 , m-M) = (359?15, -59?60, 17.01) 
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Fig. 8.— Possible tidal tails around DES J2356—5935 (Tucana III). Top left'. Spatial distribution 
of stellar objects with g < 23 mag that pass the same isochrone filter as applied to visualize 
DES J2356—5935 (Eigure 7). Red boxes highlight the positions of the tidal tails. Top right'. Spatial 
distribution of stellar objects with g — r < 1 mag and g < 23 mag not passing the isochrone filter. 
Bottom left: Average density prohles transverse to the tidal tail axis: stars passing the isochrone 
filter (red), other unfiltered stars (gray), and galaxies passing the isochrone filter (black). Eor the 
density calculation, the range of included angles parallel to the tidal tails is ±2?0, excluding the 
region within 0?3 of DES J2356—5935. Bottom right'. Binned significance diagram representing 
the Poisson probability of detecting the observed number of stars within 0?2 of the tidal tail axis 
(excluding the region within 0?3 of DES J2356—5935) for each bin of the color-magnitude space 
given the local field density. 
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Y2 DES candidates (182 kpc) and is detected as a compact cluster of RGB stars. Both the 
BHB and red horizontal branch (RHB) appear to be well populated. The distance, physical 
size (103pc), and luminosity {My = —4.5mag) of DES J0531—2801 place it in the locus 
of Local Group dwarf galaxies, comparable to Leo IV and GVn II (Belokurov et al. 2007). 
DES J0531—2801 is isolated from the other new DES systems and is likely not associated with 
the Magellanic system. 

• DES J0002—6051 (Tucana IV, Figure 10): DES J0002—6051 has the largest angular size 
of the candidates (rh = 9ll), which at a distance of 48kpc corresponds to a physical size of 
ri /2 ~ 127pc. This large half-light radius is inconsistent with the sizes of known globular 
clusters (Harris 1996, 2010 edition), thus making it very likely that DES J0002—6051 is a 
dwarf galaxy. The measured ellipticity of DES J0002—6051, e = 0.4, is also consistent with 
a galactic classification. DES J0002—6051 is found to be 27 kpc from the LMG and 18 kpc 
from the SMG. DES J0002—6051 is one of the proposed members of the Tucana group, with 
a centroid separation of 7 kpc. Measurements of the radial velocity and proper motion of 
DES J0002—6051 will provide strong clues as to whether it was accreted as part of a system 
of satellites. Similarly to DES J2356—5935 and DES J2204—4626, the MSTO of DES J0002 
—6051 is well-populated and clearly visible. Several possible member stars can also be seen 
along the HB. 

• DES J0345—6026 (Reticulum III, Figure 11): DES J0345—6026 appears similar to DES J2204 
—4626 in its structural properties, but is more distant (heliocentric distance of 92 kpc). Like 
DES J2204—4626, DES J0345—6026 can be tentatively classified as a dwarf galaxy based on 
its physical size and low surface brightness. There is some indication of asphericity for this 
object; however, the ellipticity is not significantly constrained by the DES data. DES J0345 
—6026 has a sparsely populated RGB, a few possible RHB members, and two possible BHB 
members. The MSTO for DES J0345—6026 is slightly fainter than our fitting threshold of 
g < 23 mag, thus its age is poorly constrained. 

• DES J2337—6316 (Tucana V, Figure 12): DES J2337—6316 is the second new system in the 
Tucana group and the closest to the group’s centroid (3kpc). At a heliocentric distance of 
55 kpc, DES J2337—6316 is also located 29kpc from the LMG and 14kpc from the SMC. The 
physical size (17pc) and luminosity {My ~ —1.6 mag) place DES J2337—6316 in a region of 
the size-luminosity plane close to Segue 1, Willman 1, and Kim 2. Segue 1 and Willman 1 
have sizes similar to the most extended globular clusters, but are approximately an order 
of magnitude less luminous. On the other hand, deep imaging of Kim 2 has led Kim et al. 
(2015b) to conclude that it is very likely a star cluster. DES J2337—6316 is one of the few 
new objects for which the DES data places a constraint on ellipticity. The high ellipticity, 
e = 0.7, supports a galactic classification for DES J2337—6316. Most of the high probability 
member stars for DES J2337—6316 are on the main sequence and fainter than g ~ 22.5 mag. 
As with other similarly faint objects, a few stars situated close to the lower RGB are also 
likely members. A system similar to DES J2337—6316 would have been difficult to detect in 
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DES J0531-2801 (Columba I); (asooo, ^^ 2000 , m-M) = (82?86, -28?03, 21.30) 
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Fig. 9. 


Analogous to Figure 6 but for DES J0531—2801 (Columba I). 
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DES J0002-6051 (Tucana IV); (asooo, ^^ 2000 , m-M) = (0?73, -60?85, 18.41) 
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Fig. 10 


Analogous to Figure 6 but for DES J0002—6051 (Tucana IV). 
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DES J0345-6026 (Reticulum III); (aaooo, ^^ 2000 , m - M) = (56?36, -60?45, 19.81) 
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Fig. 11. 


Analogous to Figure 6 but for DES J0345—6026 (Reticulum III). 
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SDSS with a threshold at r < 22 . 

• DES J2038—4609 (Indus II, Figure 13); DES J2038—4609 is the first of the two lower con¬ 
fidence candidates. The detection of this object comes predominantly from a tight clump 
of BHB stars at 5 ~ 22 mag and 5 — r ~ Omag. Three of the potential HB members are 
clustered within a spatial region of radius ~ IQ". Several additional HB and RGB stars are 
assigned non-zero membership probabilities by the likelihood fit, enlarging the best-fit size of 
this system. DES J2038—4609 is located at a boundary between the Y1 and Y2 imaging, and 
the effect of the deeper Y2 source detection threshold is visible when extending to magnitudes 
of S' > 23 mag. The best-fit physical size of DES J2038—4609 (181 pc) would place it among 
the population of Local Group galaxies. Deeper imaging is needed to better characterize this 
candidate since the MSTO is fainter than the current DES detection threshold. 

• DES J0117—1725 (Cetus II, Figure 14): DES J0117—1725 is the second lower confidence 
candidate, being the least luminous {My = 0 . 0 ± 0 . 68 mag) and most compact (ri /2 = 17 pc) of 
the new stellar systems. DES J0117—1725 is nearly two orders of magnitude fainter than glob¬ 
ular clusters with comparable half-light radii. At a heliocentric distance of 30 kpc, DES J0117 

— 1725 is the second nearest of the new systems. DES J0117—1725 is detected predominantly 
by its main sequence stars and has one potential HB member. The lack of strong features 
in the color-magnitude distribution of stars associated with DES J0117—1725 leads to a large 
degeneracy between the distance, age, and metallicity of this system. If determined to be 
a dwarf galaxy, DES JOl 17—1725 would be the least luminous known galaxy; however, the 
current classification of this object is ambiguous. 

There are several gaps in the Y2Q1 coverage ~ 0?1 from DES J0117—1725; however, DES J0117 

— 1725 is sufficiently compact that few member stars are expected in the region of missing 
coverage. These gaps are incorporated into the likelihood analysis and should have a min¬ 
imal impact on the fit. The region around DES J0117—1725 was imaged in only a single 
r-band exposure; however, there are two to three overlapping exposures in each of the g- and 
i-bands. We have verified that the stellar overdensity is also apparent in the g — i filter com¬ 
bination. Nonetheless, the properties of DES J0117—1725 should be interpreted with caution 
until additional imaging is acquired. 


6.2. Satellite Detection Completeness 

Many studies of the Milky Way satellite galaxy population require knowledge of the detection 
efficiency for satellites as a function of luminosity, heliocentric distance, physical size, and sky 
position in order to account for observational selection effects. To quantify our search sensitivity 
with the Y2Q1 data set, we follow the approach that Bechtol et al. (2015) applied to the DES YlAl 
coadd data. Briefly, we generated many realizations of satellite galaxies with different structural 
properties that are spread uniformly throughout the survey footprint, excluding regions of high 
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DES J2337-6316 (Tucana V); ( 02000 , ^^ 2000 , m - M) = (354?35, -63?27, 18.71) 
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Fig. 12. 


Analogous to Figure 6 but for DES J2337—6316 (Tucana V). 
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DES J2038-4609 (Indus II); (azooo, ^ 2000 , M) = (309?72, -46?16, 21.65) 
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Fig. 13.— Analogous to Figure 6 but for DES J2038—4609 (Indus II). Note that the magnitude 
threshold for the smoothed counts maps and radial profiles has been raised to g < 23 mag due to 
non-uniform imaging depth in this region. DES J2038—4609 is detected primarily by a cluster of 
potential BHB members. 
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DES J0117-1725 (Cetus II); (azooo, <^ 2000 , m - M) = (19?47, -17?42, 17.38) 



Fig. 14.— Analogous to Figure 6 but for DES J0117—1725 (Cetus II). Note that the magnitude 
threshold for the smoothed counts maps and radial profiles has been raised to g' < 23 mag due to 
incomplete coverage in this region. 
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stellar density near to the LMC. We then applied the map-based search algorithm to estimate the 
likelihood of a high-confidence detection. In the present implementation, the detection significance 
was evaluated as the Poisson probability of finding the “observed” number of stars within a circular 
region matched to the half-light radius of the simulated satellite given the effective local field density 
after the isochrone selection (see Section 3.1). Realizations that yield at least 10 detectable stars 
and produce a > 5.5(T overdensity within the half-light radius are considered “detected”. 

An essential input to the completeness calculation is the magnitude threshold for individual 
stars. For the searches described in Section 3, we used deeper magnitude thresholds (e.g., g < 
24 mag) in order to maximize the discovery potential. A more conservative threshold, g < 23 mag, 
is selected to estimate the completeness in detecting new stellar systems based on the following 
considerations. First, the Y2Q1 stellar sample is expected to be nearly complete for both Y1 and 
Y2 exposures at magnitudes g < 23 mag. Figure 3 shows that the measured stellar density in this 
magnitude range varies smoothly over the Y2Q1 footprint. At fainter magnitudes, discontinuities 
in the field density can occur in regions where exposures from Y1 and Y2 overlap (typically, more 
faint objects are detected in the Y2 exposures). Second, the field population with g > 23mag in 
high Galactic latitude regions is likely dominated by misclassified galaxies rather than halo stars, 
and therefore the assumption of isotropy on arcminute angular scales is compromised due to galaxy 
clustering. The adopted magnitude threshold matches that used in Bechtol et al. (2015), and in 
this case, the sensitivity estimates obtained in the prior work can be approximately scaled to the 
enlarged effective solid angle of the Y2Q1 data set, ~ 5000 deg^. 

In Table 4, we list the expected detection efficiencies for the DES ultra-faint galaxy candidates 
using magnitude thresholds of < 23 mag and g <22 mag, which roughly correspond to the stellar 
completeness limits of Y2Q1 and SDSS, respectively. For the brighter threshold, several of the DES 
candidates have substantially reduced detection probabilities. 

We note that all of the faint stellar systems that were discovered in the YlAl data set are also 
significantly detected using the Y2Q1 data set, suggesting that our present search has sensitivity 
comparable to previous studies. 


6.3. Total Number and Distribution of Milky Way Satellite Galaxies 

The detection of seventeen candidate ultra-faint galaxies in the first two years of DES data is 
consistent with the range of predictions based on the standard cosmological model (Tollerud et al. 
2008; Hargis et al. 2014). These predictions model the spatial distribution of luminous satellites 
from the locations of dark matter subhalos in cosmological A^-body simulations of Local Group 
analogs. Significant anisotropy in the distribution of satellites is expected (Tollerud et al. 2008; 
Hargis et al. 2014), and likely depends on the specific accretion history of the Milky Way (Deason 
et al. 2015). 

The proximity of the DES footprint to the Magellanic Glouds is noteworthy when considering 
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Table 4. Expected detection efficiencies for ultra-faint galaxy candidates. 


Name 

Mv 

Distance 


Efficiency 

Efficiency 



(kpc) 

(arcmin) 

(g < 23 mag) 

(g < 22 mag) 

DES J2204-4626 (GruII) 

-3.9 

53 

6.0 

1.00 

0.39 

DES J2356-5935 (TucIII) 

-2.4 

25 

6.0 

1.00 

0.94 

DES J0531-2801 (Coll) 

-4.5 

182 

1.9 

0.95 

0.86 

DES J0002-6051 (TucIV) 

-3.5 

48 

9.1 

0.98 

0.05 

DES J0345-6026 (Ret III) 

-3.3 

92 

2.4 

0.33 

0.05 

DES J2337-6316 (TucV) 

-1.6 

55 

1.0 

0.94 

0.01 

DES J2038-4609 (Indll) 

-4.3 

214 

2.9 

0.26 

0.00 

DES J0117-1725 (Get II) 

0.0 

30 

1.9 

0.38 

0.01 

Ret II 

-3.6 

32 

6.0 

1.00 

1.00 

DES J0344.3-4331 (Erill) 

-7.4 

330 

1.8 

1.00 

1.00 

DES J2251.2-5836 (TucII) 

-3.9 

58 

7.2 

1.00 

0.07 

Horl 

-3.5 

87 

2.4 

0.66 

0.14 

DES J0443.8-5017 (Pic I) 

-3.7 

126 

1.2 

0.85 

0.47 

DES J0222.7-5217 (Erilll) 

-2.4 

95 

0.4 

0.41 

0.04 

DES J2339.9-5424 (Phell) 

-3.7 

95 

1.2 

0.99 

0.77 

Grul 

-3.4 

120 

2.0 

0.21 

0.04 

Horll 

-2.6 

78 

2.1 

0.26 

0.01 


Note. — Expected detection efficiencies are provided for two magnitude thresholds roughly 
corresponding to the stellar completeness depths of DES Y2Q1 (g < 23 mag) and SDSS {g < 
22 mag). The detection efficiencies are estimated from multiple realizations of each candidate 
within the DES footprint using the best-fit luminosity Myi heliocentric distance, and azimuthally 
averaged half-light radius which were then analyzed with the map-based detection algorithm 
described in Section 3.1. The average detection probability ratio between the two threshold choices 
is ~ 50%, implying that roughly half of the DES candidates would have been detected if they were 
located in the SDSS footprint. 



Fig. 15.— The cumulative distribution functions of declination values for the seventeen observed 
satellites in the DES footprint (black), a uniform distribution given the shape of DES footprint 
(blue), and the best-fit model for a uniform distribution and Magellanic Cloud components (red 
dashed). See Section 6.3. 
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DES DES DES + SDSS 



Fig. 16.— Maximum-likelihood fit to the spatial distribution of satellites for a model including 
isotropic, LMC, and SMC components (Section 6.3). The horizontal axis represents the total 
nnmber of satellites detectable at DES Y2Q1 depth integrated over the entire sky (not including 
the classical dwarf galaxies and detection inefficiency near the Galactic plane). The vertical axis 
represents the fraction of these satellites associated with the Magellanic Clouds. Left: The la, 
2a, and 3 cj likelihood contours when considering only DES data. Center: la contours for different 
valnes of the LMC and SMC truncation radius (r^) and the slope of the radial prohle (a). Right: 
Contours of the likelihood function when SDSS observations are included. Red contours {la, 2a, 
3(t) assume that 50% of the satellites discovered in DES would have been detected if they were 
located in the SDSS footprint. This value matches the estimated detection efficiency from Table 4. 
We show the sensitivity of our results to the SDSS/DES detection ratio with the dark blue contour 
(80% detection ratio) and the light bine contonr (20% detection ratio). 
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an anisotropic distribution of satellites. The possibility that some Milky Way satellites are or 
were associated with the Magellanic system has been discussed for some time (Lynden-Bell 1976; 
D’Onghia & Lake 2008; Sales et al. 2011; Nichols et al. 2011), and the recent DES results have 
renewed interest in this topic (Koposov et al. 2015a; Deason et al. 2015; Wheeler et al. 2015; Yozin & 
Bekki 2015). The existence of structures on mass scales ranging from galaxy clusters to ultra-faint 
satellite galaxies is a generic prediction of hierarchical galaxy formation in the cold dark matter 
paradigm. However, specific predictions for the number of observable satellites of the LMC and 
SMC depend on the still unknown efficiency of galaxy formation in low-mass subhalos. Identifying 
satellites of satellites within the Local Group would support the standard cosmological model, and 
may provide a new opportunity to test the impact of environment on the formation of the least 
luminous galaxies (Wetzel et al. 2015). 

Based on a study of Milky Way-LMC analogs in iV-body simulations, Deason et al. (2015) 
propose a method to evaluate the probability that a given satellite was at one time associated 
with the LMC based upon the present three-dimensional separation between the two (Table 3). 
This prescription predicts that two to four of the Y2 DES candidates were associated with the 
LMC before infall, in addition to the two to four potentially associated satellites found in Y1 DES 
data. Deason et al. further suggest that a grouping of satellites near to the LMC in both physical 
separation and velocity would imply that the Magellanic system was recently accreted onto the 
Milky Way since members of this group would disperse within a few Gyr. 

As a population, the new DES satellite galaxy candidates do appear to be unevenly distributed 
within the survey footprint; fifteen are located at (52ooo < —40°, close to the Magellanic Clouds 
(Figure 3). A Kolmogorov-Smirnov test, shown in Figure 15, significantly rejects the hypothesis of 
a uniform distribution {p = 4.3 x 10“^). One possible explanation for non-uniformity is that some 
of the new systems are satellites of the LMC or SMC. Below we introduce a simple model to test 
this hypothesis. 

We began by modeling the probability distribution for detecting a satellite galaxy as a func¬ 
tion of sky position. In this analysis, we excluded classical Milky Way satellites and considered 
only satellites that have similar characteristics to those discovered by DES (i.e., are detectable at 
the DES Y2Q1 depth). The underlying luminosity function and the completeness as a function 
of luminosity, heliocentric distance, and physical size were all simplified into a two-dimensional 
probability distribution. 

We modeled this two-dimensional probability distribution with three components; (i) an 
isotropic distribution of satellites associated with the main halo of the Milky Way, (ii) a popu¬ 
lation of satellites spatially associated with the LMC halo, and (in) a population of satellites spa¬ 
tially associated with the SMC halo. The first component was uniform on the sky, and the latter 
two components were spatially concentrated around the Magellanic Clouds. We assumed spheri¬ 
cally symmetric three-dimensional distributions of satellites centered around each of the Magellanic 
Clouds having a power-law radial profiles with slope a = — 1 and truncation radius rt = 35kpc. 
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These distributions were then projected onto the sky to create two-dimensional predicted density 
maps. 

Several underlying assumptions went into this simple model. First, we assumed that the 
Magellanic components were concentrated enough that variations in the completeness with respect 
to heliocentric distance would not substantially alter the projected density maps. Second, we 
assumed that the truncation radii of the LMC and SMC are the same, but later test the sensitivity 
of our results to this choice. Third, we assumed that smaller satellites of the Milky Way do not 
have their own associated satellites. We note that a power-law radial profile with slope a = —1 can 
be viewed as a very extended Navarro-Frenk-White profile (Navarro et al. 1997). Since the DES 
footprint does not cover the central regions of the Magellanic Clouds, this shallow profile yields a 
conservative estimate on the fraction of satellites associated with the Magellanic system. 

We performed an unbinned maximum-likelihood fit of our three-component model to the dis¬ 
tribution of observed satellites. The parameters of the likelihood function are the total number 
of satellites and the relative normalizations of the three components. The likelihood is defined 
as the product of the probabilities to detect each observed satellite given its sky position. When 
presenting our results, we marginalized over the ratio of LMC component to SMC component. This 
ratio is not strongly constrained by our analysis; however, there is a slight preference for a larger 
SMC component. 

The left panel of Figure 16 shows the constraints using the sky positions of the seventeen 
DES satellite galaxy candidates. We find that the DES data alone reject a uniform distribution of 
satellites at > 3 cj confidence when compared to our best-fit model. This again implies that there is 
either a clear overdensity of satellites around the Magellanic Clouds, or that the initial assumption 
of isotropy around the main halo is incorrect. In Figure 15, we show the distribution of declinations 
for our best-fit three-component model and find good agreement with the observations {p = 0.7). 

We next tested the sensitivity of our result to the assumed distribution of satellites around 
the LMC and SMC. The central panel of Figure 16 shows the Icr contours for three different values 
of rt = {25, 35,45 kpc}, and two values of the slope, a = {-1,-2}. Despite slight changes in the 
likelihood contours, in each case > 10% of Milky Way satellites are likely to be spatially associated 
with the Magellanic Clouds. 

We further examined the satellite population by simultaneously considering the sixteen ultra¬ 
faint satellites observed in the SDSS DRIO footprint.To fully combine DES and SDSS obser¬ 
vations we would need to model the detection efficiency of SDSS relative to DES as a function of 
satellite distance, size, and luminosity. However, since the SDSS DRIO footprint only overlaps with 
the isotropic component in our model, we can again fold the complications of detection efficiency 
into a constant detection ratio. Our fiducial calculation assumes that 50% of the satellites discov¬ 
ered in DES would be detectable by SDSS, which agrees with the relative detection efficiency of 


include Pegasus III (Kim et al. 2015a) in the list of ultra-faint satellites from SDSS. 
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the seventeen DES satellites at the depth of SDSS (Table 4). The right panel of Figure 16 shows 
that the DES+SDSS constraints are in agreement with those from DES data alone, though the 
fraction of LMC and SMC satellites is more tightly constrained by including information from a 
large region widely separated from the Magellanic Clouds. The combined results imply that there 
are ~ 100 satellites over the full sky that are detectable at DES Y2Q1 depth, and that 20-30% 
of these might be associated with the Magellanic Clouds. Importantly, this prediction does not 
attempt to model the full population of Milky Way satellites beyond the DES Y2Q1 sensitivity and 
ignores the diminished detection efficiency close to the Galactic plane. 

The Pan-STARRS team has recently identified three candidate Milky Way satellite galaxies 
in their 37r survey (Laevens et al. 2015a,b). If this search is complete to the depth of SDSS over 
its full area, it is notable that so few candidates have been found. This observation alone could 
imply anisotropy without requiring the influence of the Magellanic Clouds. However, much of 
the Pan-STARRS area is located at low Galactic latitudes where the elevated foreground stellar 
density and interstellar extinction may present additional challenges. Under the assumption that 
Pan-STARRS covers the full sky with <52000 > —30°, ~2000deg^ overlap with DES Y2Q1. This 
area of the DES footprint includes two new candidates, one of which has a large enough surface 
brightness to likely have been detected at SDSS depth (DES J0531—2801, see Table 4). These two 
candidates are located in a region of the sky that would be observed at a relatively high airmass 
by Pan-STARRS and may suffer from decreased detection efficiency. 

From this analysis, we conclude that the distribution of satellites around the Milky Way is 
unlikely to be isotropic, and that a plausible component of this anisotropy is a population of 
satellites associated with the Magellanic Clouds. However, several alternative explanations for 
anisotropy in the Milky Way satellite distribution exist. For example. Milky Way satellites could 
be preferentially located along a three-dimensional planar structure, as has been suggested by 
many authors starting with Lynden-Bell (1976). This proposed planar structure encompasses the 
Magellanic Clouds and many of the classical and SDSS satellites. Pawlowski et al. (2015) suggest 
that the satellites discovered in Y1 DES data are also well aligned with this polar structure. We 
note that the Y2 discoveries presented here include several objects near the SMC and may reduce 
the fraction of objects in close proximity to the proposed plane. An additional possibility is that 
the satellites are associated with the orbit of the Magellanic System and are not isotropically 
distributed around the Magellanic Clouds themselves (Yozin &: Bekki 2015). The DES footprint 
covers only a fraction of the region surrounding the Magellanic Clouds and additional sky coverage 
may yield more satellites with similar proximity to the Magellanic system and/or help to distinguish 
between these various scenarios. Measurements of the relative motions of the satellites and further 
theoretical work will also help clarify the physical relationships between these stellar systems. 
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7. Conclusions 

We report the discovery of six new ultra-faint galaxy candidates in a combined data set from the 
hrst two years of DES covering ~5000deg^ of the south Galactic cap. Two additional candidates 
are identihed in regions with incomplete or non-uniform coverage and should be viewed with lower 
confidence until additional imaging is obtained. The new satellites are faint (My > —4.7 mag) and 
span a wide range of physical sizes (17pc < r^i 2 < 181 pc) and heliocentric distances (25kpc < 
Hq < 214 kpc). All are low surface brightness systems similar to the known ultra-faint satellite 
galaxies of the Milky Way, and most possess physical sizes that are large enough (> 40 pc) to be 
provisionally classified as galaxies. Spectroscopic observations are needed to better understand 
and unambiguously classify the new stellar systems. A total of seventeen confirmed and candidate 
ultra-faint galaxies have been found in the first two years of DES. Roughly half of the DES systems 
are sufficiently distant and/or faint to have eluded detection at survey depths comparable to SDSS. 

The DES satellites are concentrated in the southern half of the survey footprint in proximity 
to the Magellanic Clouds. In addition, we hnd three satellites clustered in a Tucana group, each of 
which is within < 10 kpc of the group centroid. We hnd that the DES data alone exclude {p < 10“^) 
an isotropic distribution of satellites within the Milky Way halo, and that the observed distribution 
can be well, although not uniquely, explained by a model in which several of the observed DES 
satellites are associated with the Magellanic system. Under the assumption that the total satellite 
population can be modeled by isotropic distributions around the Milky Way, LMC, and SMC, 
we estimate that a total of ~ 100 ultra-faint satellites with comparable physical characteristics to 
those detected by DES might exist over the full sky, with 20-30% of these systems being spatially 
associated with the Magellanic Clouds. 

Milky Way satellite galaxies are considered a unique population for studying the particle 
nature of dark matter due to their proximity, characteristically large mass-to-light ratios, and lack 
of intrinsic astrophysical backgrounds. Although the dark matter content of the new stellar systems 
has not yet been spectroscopically confirmed, the possibility that some are dark matter dominated 
galaxies makes them interesting targets for indirect dark matter searches via gamma rays. No 
gamma-ray sources from the third Fermi-LAT source catalog (3FGL; Acero et al. 2015) are located 
within 0?5 of any of the new candidates. A follow-up gamma-ray analysis of these ultra-faint galaxy 
candidates will be presented separately. 

Future seasons of DES will not significantly increase the DES sky coverage. However, there 
will be considerable gains in depth from the expected ten tilings per hlter relative to the two to 
four in the Y2Q1 data set. This increased depth will expand the effective volume of the survey. 
In addition, we expect the survey uniformity to improve with the coming seasons, resulting in a 
cleaner list of seeds with fewer false positives. Star-galaxy separation may become a limiting factor 
at the future depth of DES, motivating the development of improved classification algorithms (e.g., 
Fadely et al. 2012; Soumagnac et al. 2015), as well as alternative search strategies, involving, for 
instance, the time domain (Baker Sz Willman 2015) or stellar velocities (Antoja et al. 2015). While 
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it is likely that the most conspicuous Milky Way satellites in the DES footprint have been found, 
the most exciting phase of Milky Way science with DES is likely still to come. The wide area and 
growing sensitivity of DES will enable the discovery of dwarf galaxies that are fainter, farther, and 
lower surface brightness. 
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